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ABSTRACT 

We study the self-similar collapse of an isothermal magnetized rotating cloud 
in the ideal magnetohydrodynamic (MHD) regime. In the limit of small distance 
from the accreting protostar we find an analytic solution that corresponds to 
free-fall onto a central mass point. The density distribution is not spherically 
symmetric but depends on the mass loading of magnetic field lines, which can be 
obtained by matching our inner solution to an outer collapse solution previously 
computed by Allen, Shu & Li. The concentration of magnetic field trapped by the 
central mass point under field-freezing, independent on the details of the starting 
state, creates a split monopole configuration where the magnetic field strength 
increases as the inverse square of the distance from the center. Under such 
conditions, the inflow eventually becomes subalfvenic and the outward transfer 



-2- 



of angular momentum by magnetic braking very efficient, thus preventing the 
formation of a centrifugally supported disk. Instead, the azimuthal velocity of 
the infalhng gas decreases to zero at the center, and the gas spirals into the star. 
Therefore, the dissipation of dynamically important levels of magnetic field is 
a fundamental requisite for the formation of protoplanetary disks around young 
stars. 

Subject headings: lSM:clouds — ISM: magnetic fields — magnetohydrodynamics 
— planetary systems: protoplanetary disks — stars: formation 



1. Introduction 

Since the realization that magnetic fields exist in interstellar clouds, there has been 
the potential for a severe magnetic-fiux problem in star formation. If field freezing were to 
apply (the ideal MHD conditions), newly formed solar- type stars should have surface fields 
of ~ 10^ G (e.g., Mestel & Spitzer 1956), which is several thousand times larger than given 
by the observations of T Tauri stars (e.g., Johns-Krull, Valenti, & Saar 2004). Magnetic-field 
dissipation is crucial to the resolution of this contradiction. In this paper, we demonstrate 
that it is also needed to avoid an extremely efficient magnetic braking that would prevent 
the formation of the observed protoplanetary disks around young stars. 

Allen, Shu & Li (2003; hereafter A03a) performed informative numerical simulations of 
the collapse of isothermal magnetized toroids, which are axially symmetric equilibrium con- 
figurations with spatially uniform mass-to-fiux ratio (such configuration is called isopedic). 
These are pivotal states at time t = 0, in the sense that they represent the idealized state of 
a molecular cloud core at the instant of formation of a central density cusp, and separate the 
pre-protostellar phase of evolution from the protostellar accretion phase (Li & Shu 1996). 
The structure of these toroids is completely determined by the single parameter Hq, that 
is the fraction of support provided against self-gravity by poloidal magnetic fields relative 
to that of gas pressure (Li & Shu 1996). The collapse of singular isothermal toroids occurs 
inside-out in a self-similar manner (A03a), in which matter collapses toward the central point 
mass dragging magnetic field lines along with it. 

Later, Allen, Li & Shu (2003; hereafter A03b) included angular momentum character- 
ized by a fiat rotation curve in these models. The formation of a centrifugally supported 
disk around the central star did not occur for any reasonable nonzero value of the magnetic 
support parameter Hq. They concluded that, in ideal MHD conditions (i.e., no magnetic 
field dissipation), the efficiency of magnetic braking was so high as to prevent disk forma- 
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tion. Simulations with limited numerical resolution can be suggestive, but proof requires an 
analytic or semi-analytic demonstration. Analytic or semi-analytic results have the added 
advantage that they can be the starting points for additional development, such as providing 
a background equatorial collapse against which a bipolar outflow might simultaneously be 
taking place along the poles (Lizano et al. 2006, in preparation), or for an extension into the 
non-ideal regime (Shu et al. 2006). 

Motivated by these considerations, in this paper we determine analytically the asymp- 
totic behavior of the density, velocity and magnetic field in the self-similar collapse of an 
isothermal magnetized rotating cloud. We find that at small scales the density and poloidal 
velocity field in the collapse solution approaches free-fall on a split magnetic monopole con- 
figuration, with the azimuthal velocity decreasing to zero at the origin. This result proves 
that, without magnetic field dissipation, a centrifugally supported disk around the accreting 
central star cannot form. 

The organization of the paper is as follows. In §2 we discuss the relevant ideal MHD 
equations. In §3 we put these equations in non-dimensional self-similar form. In §4 we write 
the self-similar equations in a streamline formulation, in terms of the Grad-Shafranov and 
Bernoulli's equation, and take the limit for small values of the self-similar variable x = r/at, 
where r is the spherical radius, a is the sound speed of the cloud and t is the time variable. In 
§5 we show some numerical examples of our inner solution for the case of no rotation, using 
the numerical simulations of A03a to obtain an approximation for the mass loading of field 
lines. In §6 we consider the effects of rotation and show that the azimuthal velocity vanishes 
at the origin, i.e., the efficient magnetic braking prevents the formation of a centrifugally 
disk around the star. Finally, in §8, we discuss the implications of our work and summarize 
our conclusions. 



2. Ideal MHD equations 

The equations governing the collapse of an isothermal magnetized cloud are the equation 
of continuity 

|^ + V-(pu) = 0, (1) 

where p is the density and u is the velocity of the gas; the equation of motion 

^ + (V X u) X u + -VluP = -— Vp -W + - — (V X B) X B, (2) 

dt ^ ' 2 ' ' p Anp^ ' ^ ' 



-4- 



where V is the gravitational potential and B is the magnetic field, the induction equation 
in the ideal MHD limit 

^ = V X (u X B), (3) 

the condition of no monopoles, 

V • B = 0, (4) 

and Poisson's equation 

V V = AttGp, (5) 

where G is the gravitational constant. 

The geometry of the collapse is shown in Figure 1. A separatrix indicated in dashed 
lines separates field lines which have been pulled into the origin by the infiow from field lines 
tied to matter in the infalling envelope (see, e.g.. Figures 4, 6 and 7 of A03a). We will study 
only the region inside the separatrix surface that is close to the origin (within the dotted 
circle in Figure 1). 



3. Self-similar variables 

Since the singular isothermal toroids defining the pivotal t = states are self-similar, 
we look for collapse solutions depending only on the similarity coordinate x, defined by 

r = atx. (6) 

We introduce appropriate nondimensional (reduced) variables defined as 

P(r,t) = ^a(x), (7) 

u(r, t) = av(x), (8) 

V{r,t) = a'^{^), (9) 

B(r,t) = ^b(x). (10) 



Inserting the definitions (6)-(10) into eq. (l)-(3), we obtain the set of nondimensional 
equations governing the self-similar collapse, 

(v-x) ■ V« + aV- V = 2a, (11) 

11 1 

-(x ■ V)v + (V X v) X V + -V|vp = Va - V^ + -(V x b) x b, (12) 

2 a a 

-(x- V)b = b + V X (vx b), (13) 

V-b = 0, (14) 

V^^ = a. (15) 
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4. The inner limit 

In the limit |x| ^ 1 we approximate for any vector or scalar variable f, 

f ~(x-V)f <(v-V)f. (16) 

With these approximations, eq. (11)-(13) become 

V ■ (av) = 0, (17) 

11 1 

(V X v) X V + -V|v|^ = Va - V^ + -(V x b) x b, (18) 

2 a a 

V X (v X b) = 0. (19) 

These equations for the self-similar variables in the inner region are formally equivalent to 
the dimensional equations of steady, ideal MHD flow. The time-dependent self-similar flow 
becomes quasi-steady when the time it takes the flow to cross the inner region becomes small 
compared to the evolutionary time governing changes of the infall envelope that feeds mat- 
ter and fleld into the inner regions. The latter condition holds independent of the starting 
assumption of self-similarity of the overall collapse (see e.g. Galli & Shu 1993b), and there- 
fore our results have greater generality than its formal derivation here. In the steady-flow 
approximation, we may profltably introduce a streamline formalism (e.g, Shafranov 1966; 
Shu al. 1994). 

In the following we adopt spherical coordinates {x,6,ip), and assume that the system 
is axisymmetric {d/dip = 0). Then eq. (17) is satisfled introducing a stream function ip{'x.) 
such that 



Vx 



xsmu 



(20) 



where the components of the poloidal velocity Vp are 



and ve = : — ;^t— • (21) 



ax"^ sin 6 06 ax sin 6 dx 

Reflection symmetry requires ip to be an odd function of 6 with respect to ^ = 7r/2, therefore 



ij (x, I) = 0. (22) 



Imposing that the integral over solid angle of the mass flow —r'^pUr gives the accretion rate 
M = rrio^l + Ho)a^/G, where mo = 0.975 is the reduced mass at the origin (Li & Shu 1997), 
we obtain the condition 

ij{x,(}) = mo{l + Ho). (23) 
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Thus, the value of iIj/iIj{x, 0) is the fraction of the accretion rate carried by the streamhnes 

from —ip to ip. 

Eq. (19) imphes that v || b. We define the function f3, related to the mass loading of 
field lines, by 

b = /3av. (24) 

Using this definition, the condition of no monopoles, eq. (14), and the equation of continuity, 
eq. (17), one obtains 

V ■ V/9 = 0, (25) 

i.e., [3 = f3{ip), is constant along streamlines. We note that (3 is antisymmetric with respect 
to the equatorial plane (by convention, f3 negative corresponds to the upper hemisphere). 

If we dot the force equation (18) with v we obtain Bernoulli's equation 

-\^r\^ + lna + ^ = H{■^lJ), (26) 

where -&('?/') is Bernoulli's function which is constant on streamlines. Substituting this equa- 
tion and using the definition of /9, eq. (24), in the equation of motion, eq. (18), we obtain 
after collecting terms, 

E'\/i) + [V X (1 - a/3=^)v + V/3 X /5av] x v = 0, (27) 

where a prime denotes differentiation with respect to ■?/'. We now decompose this equation 
along &^p and V-?/'. 

After some algebra, one can write the poloidal component of eq. (27) as 

a' + — ^ - ^J' - o^m-^f) V^ = 0, (28) 

ax sm V X sm U / 

where 

cje^ = V X (1 - a/32)vp (29) 

is the total vorticity. For nontrivial solutions of this equation to exist, we require 

H' + — ^ - ^3' - /3/5'a|vp = 0, (30) 

ax sm 6 x sm 6 

which is a form of the Grad-Shafranov (G-S) equation that expresses the condition of force 
balance across field lines. We can write the total vorticity as 



1 

UJ = 



xsin^ 






(31) 



-7- 



where S is the Stokes operator 

d'^ip 1 d'^ip cot 6 dip 

and A is the Alfven discriminant 

A='-^^. (33) 

a 

Then, the G-S equation, (30), can be written as 

AS(ij) = aH'x^ sin^ e + irVa -Vij - \ jj' + pp'lVipl^. (34) 

a^ A 

The toroidal component of equation (27) is 

V • (Vj + af3v^Vf3) = 0, (35) 

where we have defined the total (gas plus magnetic field) specific angular momentum 

j = xsm9(l — aP'^)v^. (36) 

Since f3 is constant along streamlines by eq. (25), eq. (35) reduces to 

V ■ Vj = 0, (37) 

stating that j = j{iIj) is constant along streamlines. Thus, the toroidal velocity is given by 

''~ xsinOil-aP^)- ^ ' 

Notice the formal singular point that arises when af3'^ = 1, which is the condition of the 
poloidal velocity reaching the Alfven speed. In fact, the smooth crossing of this surface (with 
a continously varying Vy,) is made in a time-dependent way; as will be shown in Sect. 5, the 
poloidal velocity is distinctly subalfvenic {a (3"^ ^ 1) in the steady part of the fiow. Finally, 
Bernoulli's eq. (26) can be rewritten as 

|V^|^+(^) +2aVsin2^(lna + ^-iJ) = 0. (39) 

Thus, in the inner limit, the ideal MHD equations of the problem (17)-(19), can be 
reduced to the G-S and Bernoulli's equations (eq. 34 and 39, respectively). These need to be 
solved for the stream function ip and the density a (or the Alfven discriminant A), given the 
functions H{iIj), j{ip), and i3{iIj), which are quantities conserved along streamlines. Then, 
the toroidal components of the magnetic field and the velocity can be obtained from eq. (24) 
and (38), respectively. 



a{x,e) = 


'moil+HoY 
2 


L/2 


and 


ij{x,e)=mo{l + Ho)f{e), 


where Q{9) is normalized so that 


/' 

Jo 


r/2 

Q{e)smede = i. 



5. Case of no rotation (j = 0) 

We first consider the case of no rotation, and we look for a separable solution in spherical 
coordinates in which all variables have a power-law dependence on x. One can show that 
the density and the stream function must have forms that represent magnetically-modulated 

free-fall solutions: 

r™ n , tj Ml/2 

Q{9), (40) 

(41) 

(42) 

The angular part of the stream function, f{0), is normalized so that /(O) = 1 and /(7r/2) = 0, 
i.e., the streamlines between those falling in equatorially and along the pole carry the entire 
mass infall rate of the upper hemisphere. 

To find the unknown functions Q{6) and f{6) we need to balance the highest order 
terms in eqs. (34) and (39). The dominant terms that balance in the G-S equation, eq. (34), 
are then 

-f3'S{^) ^ /3/5'|V^p, (43) 

where the Alfven discriminant, A ^ — /3^, corresponds to highly subalfvenic poloidal flow. 

In Bernoulli's equation the dominant terms that balance are the poloidal kinetic energy 
and the gravitational potential energy. With the gravitational potential of a point mass at 
the origin, \1' = —7710(1 + Ho)/x, we get 

\Vij\'^ ^2moil + Ho) a^xsin^ 9. (44) 

Substituting the definitions (40) and (41) in this latter equation, one gets 

^ = -g(^)sin^. (45) 

The G-S equation (43)now becomes a second-order, nonlinear, ordinary differential equation 

for f{e), 
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Combining the Bernoulli equation (45) and the G-S equation (46), one obtains 

df j3 

which has the solution 

QU) = ^. (48) 

where j3 is an integration constant. Furthermore, f{d) has the parametric solution 

cos^(/) = i / /3(/)rf/, (49) 

P Jo 

where we have used the B.C. /(7r/2) = 0. Moreover, the B.C. /(O) = 1, implies that 

'P= f PU) df. (50) 

The physical meaning of the above results is simple. Equation (49) implies that the 
poloidal field lines spread so as to distribute themselves uniformly in polar angle, i.e., the 
magnetic configuration takes the shape of a split monopole, as anticipated by Galli & Shu 
(1993a; see their Figure 2) and confirmed later by Li & Shu (1997) for the collapse of an 
infinitesimally thin disk. This happens because the field lines trapped in the central mass 
point become so strong that the mass infiow, being highly subalfvenic, is unable to bend them 
from the zeroth-order vacuum field configuration. In turn, the matter is constrained by field 
freezing to flow along the radial field lines, and the density distribution is that appropriate 
for free-fall along each radial path (eq. 41), but it departs from a homogenous angular 
distribution in a manner that depends only on the initial mass-to-fiux loading (eq. 48). 

Given the above interpretation, it is useful to define the magnetic flux function 

*(r>^) = |^</>W, (51) 

from which we may obtain the poloidal fleld B: 

The components of the poloidal magnetic fleld, in nondimensional units, are given by 

and bg = ——— — . (53) 



x^ sin 9 89 X sin 9 dx 
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Comparing these definitions with equations (21), and using the definition oi (3, eq. (24), one 
sees that there is a simple relation between the fiux and the stream function given by 

d(f) = (3dilj (54) 

Thus, the parametric equation for ip, eq. (49), simply states that 

= 0^(1- cos ^), (55) 

where 0* = 0(7r/2) (just above the midplane) is the nondimensional magnetic flux trapped 
in the central source. In fact, 0^ labels the flux of the separatrix (see Figure 1). Moreover, 
the normalization eq. (50) and eq. (54) give 

Since every fleld line that lies interior to the separatrix threads the mass point at the center, 
one has the immediate identiflcation that 

A. = -r' = '"°'^ + ^°> . (57) 

where A^, is the nondimensional mass-to-flux ratio of the central star deflned as 

A.^?^^^. (58) 

$(r,7r/2) ^ ' 

To summarize, the solution for small x consists of streamlines which are radial but non- 
uniformly distributed in 9. The 9 dependence of both -0 and a is determined by the mass 
loading of fleld lines represented by i3{iIj), assumed to be known from the matching to the 
outer collapse solution. According to eqs. (21) and (45), the radial velocity is just free-fall, 
i.e., 

'2moil + Ho)Y^^ 



X 

By eq. (24), the magnetic fleld is that of a split monopole at the the origin. 



(59) 



K = % (60) 

X'^ 

as inferred by A03a,b. In physical units, the magnetic fleld of the split monopole is 

Br = 0.^, (61) 
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and is maintained by a current sheet in the equatorial plane, with current density in the 

azimuthal direction 

cB 

J, = ^5(^-7r/2). (62) 

Since the radial velocity increases as Vx oc x~^/^, while the (nondimensional) Alfven 
speed increases much faster, va = bx/cx^^'^ oc a;"^/^, the flow asymptotically becomes highly 
subalfvenic, consistent with our assumption for A in eq. (43). 



6. Case with rotation (j ^ 0). 

In the case j = 0, the balance across and along streamlines is achieved by terms of order 
x~^ (see eq. 43 and 44). Multiplying eq. (34) by A, and eq. (39) by A"^, one can show that 
the highest-order balance of the G-S and Bernoulli equations for j 7^ is unchanged. Thus, 
the non-rotating solution for the density a and the radial velocity and magnetic field, Vx 
and bx, derived in the previous Section remains valid for x <^ 1, for any value of j. This is 
consistent with the findings of A03b that the isodensity contours and poloidal field (stream) 
lines are the same to zeroth approximation in the rotating and non-rotating case for small 

X. 

Using eq. (37) together with the equation of continuity, eq. (17), the conservation of 
total angular momentum, V ■ (avj) = 0, can be written as 

V ■ {a\pjg) = V • (hpb^xsmO) (63) 

where the gas specific angular momentum is jg = xsin^f,^. This equation can be integrated 
over the upper hemisphere as 

/ Oijg^p ■ (iS = / xsinOb^hp ■ dS, (64) 

where dS = 2'kx'^ sin 6d6 ex- The left hand side is rate of change of the gas specific angular 
momentum and the right hand side is the magnetic torque. Given our inner solution, the 
integral on the left hand side is oc xv^p and the magnetic torque is oc a:~^/^f<^. Thus, the 
magnetic torque eventually dominates, decreasing the net gas angular momentum and forcing 
the azimuthal velocity to deviate from the Vy, oc x^^, behavior valid in the absence of 
magnetic torques. As discussed above, for small x the flow becomes subalfvenic {af3'^ ^1). 
Thus, eq. (38) for the azimuthal velocity, together with eq. (40), implies that along each 
streamline v^ oc x^^"^, i.e., the azimuthal velocity decreases to zero at the origin. The 
long lever arm associated with the divergently strong split monopole field yields highly 



-12- 



efficient magnetic-braking, causing the gas to spiral into the central star without forming a 
centrifugally supported disk, as suggested by A03b. 

Finally, from eq. (24), one can obtain the toroidal magnetic field component b^ oc x~^, 
that increases with decreasing radius more slowly than the poloidal component b^ (eq. (60)). 
Thus, the winding of the field {b^/b^ oc x) goes to zero as a; ^ 0, where the field geometry is 
dominated by the split monopole at the origin. Thus, in our treatment, reconnection of the 
poloidal magnetic field because of the sheet current that prevails in the magnetic midplane 
is more important than reconnection of the toroidal component that results because of the 
spinup of inflowing matter. In realistic circumstances, this state of affairs may change after 
the poloidal field is nearly entirely dissipated and a centrifugal disk is formed that can wind 
up relatively weak fields into a predominantly toroidal configuration. 

Figure 2 summarizes the results of this Section and illustrates qualitatively the behavior 
of the azimuthal velocity in the equatorial plane of a collapsing magnetized cloud. We 
stress that different choices of the mass-loading of field lines specified by the function f3{ip) 
would merely change the details of when and where the maximum rotational velocity is 
reached in the inflow solution, not the fact that the inflowing matter eventually loses all 
its specific angular momentum. In the next section, adopting a specific collapse model, 
we will show quantitatively that the region where magnetic braking by the split monopole 
becomes dominant over angular momentum conservation is generally larger than the size of 
a centrifugally supported disk formed in the absence of magnetic torques. 



7. Matching to the outer collapse solution 

We can match our inner collapse solution a; <^ 1 to the outer collapse solution via the 
mass-loading function, P{ip), and the specific total angular momentum, j{ip). We use the 
simulations of A03a to perform this asymptotic matching. The exercise is nontrivial since 
collapse simulations with very high grid resolution are required to provide the "inner limit 
of the outer solution," onto which we wish to match the "outer limit of the inner solution." 
For our purposes, there is a further problem with the boundary condition imposed by Allen 
and co-workers on the radial field, namely fe^; = at ^ = 7r/2. This boundary condition, 
appropriate to the outside of the separatrix, is in conflict with the magnetic field configuration 
inside of the separatrix (that becomes the split monopole solution at small x) depicted in 
Figure 1. In fact, unlike the analytic work, it is numerically impossible to allow a current 
sheet at the equator (where bx is discontinuous at the midplane) since this configuration is 
unstable to magnetic reconnection via numerical diffusion (see, e.g.. Figure 4 of Galli & Shu 
1993b). In spite of these drawbacks, we use the numerical [3 obtained from the simulations of 
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A03a to apply to our models, trusting the conservation principle associated with field freezing 
to yield the correct mass-to-fiux loading for all colatitudes except close to the equator. 



In Figure 3 the logarithm of |/3| is plotted as a function of the normalized fiux, ^ = 
for different values of Hq increasing from bottom to top. The points are the results of the 
simulations. The broken lines are fits given by 

-Pho = 6exp(-c^)r''" (1 - ey + Ho, (65) 

where b = Hq"^ exp{3.3Ho) and c = 1.6 + 2.3Hq. These broken lines yield a good represen- 
tation of the numerical data for all angles except those close to the equator, ,^ = 1. The 
numerical f3 goes to zero due to the boundary condition discussed above (see eq. 24). The 
analytic solution requires [3 to be finite and /3' = at the midplane; otherwise the density 
gradient will become discontinuous. 

Note that the function j3 diverges at the pole [6 = 0) because in the starting toroids the 
density there goes to zero. Finally, since the numerical simulations give /5(0), not f3{ip), we 
use the fact that dip = (f)-^d$,/f3{C,), to obtain from eq. 54) and (57) 

The numerical values of A^ can be determined only approximately, owing to the effects of 
the boundary condition at the midplane on the numerical values of f3. Roughly, we obtain 
A^ ^ 3-4 for Ho = 0.125, A^ ^ 2.4-2.7 for Hq = 0.25, A^ ^ 1.2-1.7 for Hq = 0.5, and 
A^ = 1.0-1.3 for Hq = 1. These values are closer to the mass-to-flux ratio A^. than they are 
to A as tabulated in Table 1 of Li & Shu (1996). The former applies for the initial state of 
an isopedic toroid when one integrates for the mass along flux tubes for given Hq only to a 
spherical surface of radius r, whereas the latter applies if one integrates all the way out to oo. 
For given trapped flux, we can expect mass accumulation to behave in the latter case only 
in the case of complete flattening, Hq —^ oo. For more modest values of Hq, we can expect a 
more nearly spherical reach for both mass and flux when the conflguration collapses, yielding 
the closer correspondence of A^ to A^. 

In contrast, T Tauri stars have observed values of A^ ~ 10^-10^ (see e.g. Johns-Krull 
et al. 2004), implying that the assumption of fleld- freezing must break down at some point 
in the star formation process (Shu et al. 2006). As expected, according to our ideal MHD 
solution the magnetic fleld trapped in the central protostar exceeds observed stellar flelds by 
three-four orders of magnitude. The four panels of Figure 4 show our inner collapse solutions 
in self-similar coordinates for the four values of Hq given above. The heavy solid contours 
in each panel are isodensity contours and the radial lines are the magnetic fleld lines (which 
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coincide with the streamhnes) . The arrows show the velocity field at different radii. The 
dynamic pseudodisks that result may bear comparison with observations of young protostars 
before they have developed full disks (Hogerheijde 2004). 

Quantitatively, the condition for subalfvenic radial fiow ajS"^ > 1, necessary for strong 
magnetic braking according to the results described in the previous section, is valid in a 
region x < xa{0) where x^ is the nondimensional Alfven radius. With our asymptotic 
expressions for v^ and b^ (eq. 59 and 60) and our fitting formula for P (eq. 65), it is possible 
to estimate the value of xa on the midplane {9 = 7i/2), obtaining xa ~ 0.085, 0.18, and 0.45 
for Hq = 0.125, 0.25 and 0.5. These values can be compared with the expected disk radius 
for the collapse of a non-magnetized Toomre-Hayashi toroid, Xd ~ 0.25fQ, where vq is the 
(uniform) rotation velocity in units of the isothermal sound speed a (A03b). The magnitude 
of fo for molecular cloud cores ranges from 0.03 to 0.4, with typical value vq ~ 0.1 (Goodman 
et al. 1993), in which case Xd ~ 0.0025. Thus, Xd <^ xa even for weakly magnetized clouds 
{Ho = 0.125-0.5), and the formation of a centrifugally supported disk cannot take place. 
Even for the weakly- magnetized, fast-rotating cloud model with vq = 0.5 and Hq = 0.25 
considered by A03b (see their Fig. 7), the radius of the disk that would be formed in the 
absence of magnetic torques, Xd ~ 0.063, is about 3 times smaller than the Alfven radius, 
Xa ~ 0.18, explaining why A03b did not observe the formation a centrifugal disk even in 
this extreme model. 



8. Comparison with previous work 

The process of magnetic braking has been usually considered in the context of the quasi- 
static evolution of a molecular cloud core (e.g. Gillis, Mestel & Paris 1974,1979; Mouschovias 
& Paleologou 1979). Recent semi- analytical calculations and numerical simulations have 
extended these works to follow the increase of density towards the formation of an optically 
thick core (Tomisaka 2002, Machida et al. 2005a,b) and through the formation of a point-mass 
into the protostellar accretion phase (Li & Shu 1997; Ciolek & Konigl 1998; Contopoulos, 
Ciolek & Konigl 1998; A03a,b). Krasnopolsky & Konigl (2002) extended the semi-analytical 
calculations by Contopoulos et al. (1998) to study the collapse of magnetized, rotating, thin 
disks surrounded by an external medium with specified physical properties. They found that 
the decoupling of matter and magnetic field (by ambipolar diffusion) is carried by an outward 
propagating shock, as first suggested by Li & McKee (1996), and the infalling rotating 
material is halted by a centrifugal shock inside which a centrifugal disk is established. The 
formation and radial extent of the centrifugal disk depend on two main parameters: the ratio 
6 of the azimuthal and vertical components of the magnetic field at the disk surface and the 
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Alfven speed in the external medium, both controlhng the importance of magnetic braking 
during the collapse. In the ideal MHD case (no ambipolar diffusion), for models with a large 
winding of the magnetic field {S = 10) they obtain complete braking, i.e. no centrifugal disk 
is formed. As in our case, the horizontal flow approaches free-fall onto the central mass and 
the magnetic field of the protostar is that of a split monopole at the origin. However, our 
model differs from theirs in that our collapsing cores are 3-D with axial symmetry and we do 
not include an external medium. The magnetic braking is due to the field lines that connect 
the inner region, a; <^ 1, with the rest of the collapsing cloud and thus the magnetic field 
winding and magnetic braking are calculated self-consistently (through the functions j and 
(3). In addition, we find that the radial component of the field largely dominates over the 
azimuthal component at small radii, so the winding of the field is never as severe as assumed 
by Krasnopolsky & Konigl (2002). 



9. Summary and conclusions 

We have examined the properties of the inner regions of an ideally magnetized, isother- 
mal, rotating cloud undergoing gravitational collapse. The density distribution and the 
poloidal velocity field in the collapse solution approach free-fall forms, while the trapped 
magnetic field rooted in the central protostar acquires a split monopole configuration. An 
analytic solution can be obtained for the central regions if numerical simulations are used 
to obtain the loading function that determines how much mass is coming down a given field 
line. For illustrative purposes, we have adopted the mass loading function (3 from the nu- 
merical simulations of A03a. We note, however, that the properties of our collapse solution 
discussed in Sect. 5 and Sect. 6 are general, and do not depend on these particular numerical 
simulations. 

The long lever arm associated with the strong field is able to brake the rotation of the 
infalling gas so efficiently that the azimuthal velocity goes to zero at the origin. In other 
words, the infalling gas spirals directly into the center without the formation of a centrifugally 
supported disk. Although this conclusion rigorously holds only for the particular class of 
collapse models considered in this paper, the assumption of self-similarity helps to brings out 
the dependence of our results on the fundamental physical parameters of the problem, like 
the cloud's mass-to-flux ratio and rotation rate, suggesting a broader validity of our main 
conclusions. Thus, the catastrophic magnetic braking discussed in Sect. 6, for example, does 
not depend on the particular choice of the outer collapse solution but follows directly from 
the ideal MHD equations. This work therefore proves that magnetic fleld dissipation is a 
crucial ingredient to allow the formation of the circumstellar disks observed around young 
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stars. It also yields departure points for calculations which include outflow in the presence 
of inflow (Lizano et al. 2006, in preparation) or modifications introduced by the effects of 
finite resistivity (Shu et al. 2006). 
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Fig. 1. — Geometry of the magnetic field in a collapsing magnetized cloud. The separatrix 
{dashed curve) separates field lines which have been pulled into an accreting protostar (at 
the origin of the coordinate system) from field lines tied to matter in the infalling envelope. 
A dashed circle indicates the inner region with nondimensional radius a; <^ 1 studied in this 
paper, where magnetic field lines are asymptotically radial. 
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Fig. 2. — Schematic behavior of the azimuthal velocity u^ {solid curve) in the equatorial 
plane of a collapsing cloud as function of the distance r from the central protostar. In the 
absence of magnetic torques, a centrifugally supported disk is formed inside a radius Vd, where 
the azimuthal velocity increasing as r~^ becomes equal to the keplerian velocity around the 
protostar (oc r~^'^). In contrast, the magnetic torque associated with the split monopole field 
of the central protostar reduces the angular momentum of the infalling gas, constraining the 
azimuthal velocity to decrease as r^/^ at small radii. Magnetic braking becomes dominant 
over angular momentum conservation when the infall velocity Ur becomes smaller than the 
local Alfven speed va- 
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Fig. 3. — The function f3 related to the mass loading of field lines (eq. 24), plotted as 
function of the normalized flux 0/0* = 1 — cos^. The triangles show the function P computed 
numerically by A03a for Hq = 0.125, 0.25, 0.5 and 1; the dashed curves are analytical fits 
given by eq. (65). 
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Fig. 4. — Inner collapse solutions (valid asymptotically for a; <^ 1) obtained for the values of 
the parameter Hq and the functions f3 shown in Figure 3. The horizontal and vertical axis in 
each panel are the cylindrical self-similar coordinates, w = xsin^ and z = xcos^. The heavy 
solid contours in each panel are isodensity contours for the levels a = 100, 30, 10, 3, 1 and 
0.3. The thin solid lines are the magnetic field lines (which coincide with the streamlines). 
The arrows show the velocity field at different radii, logarithmically spaced in the velocity 
interval [2.5, 6]. 



